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Abstract 

Using exact diagonalization techniques we study the dynamical response of the 
anisotropic disordered Heisenberg model for systems of S = 1/2 spins with infinite 
range random exchange interactions at temperature T = 0. The model can be 
considered as a generalization, to the quantum case, of the well known Sherrington- 
Kirkpatrick classical spin-glass model. We also compute and study the behavior 
of the Edwards Anderson order parameter and energy per spin as the anisotropy 
evolves from the Ising to the Heisenberg limits. 
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Many real materials display a spin-glass phase at low enough temperature. 
This exotic state is characterized by a frozen configuration of local magnetic 
moments following a random spatial pattern, in such a way that no net macro- 
scopic magnetization is produced. 

The manifestations of this state were first observed in the seventies in systems 
with magnetic impurities in a metallic host (examples are Cu-Mn, Ag-Mn, Au- 
Mn and Ag-Mn). There the magnetic moments experiment long range RKKY 
interactions, but spin-glass behavior was also subsequently observed in the 
insulating compound Eu^Sri.^S with competing ferromagnetic and random 
antiferromagnetic exchange interactions [1]. More recently, spin glass phases 
have also been observed in the bi-layer kagome SrCr 8 Ga40ig [2], in the py- 
rochlore structure Li I .Zni_ x V204 [3], in the dipolar magnet LiHo x Yi_ a; F4 [4] 
and in the enigmatic high Tc compounds Lai_ x Sr a ,Cu204 [5]. 
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Competing interactions and frustration in combination with randomness have 
been identified as the basic ingredients to drive a system towards a glassy 
state. While the concept of spin is purely quantum, it is usually stated that 
quantum fluctuations are not important to describe the spin glass physics. 
However, the relevance of quantum effects is beginning to be identified and 
emphasized in experimental [6,7,8,9] and theoretical work [10,11,12,13]. 

Exact diagonalization (ED) techniques proved to be very useful and efficient to 
investigate the dynamical properties of strongly correlated systems in general 
[5]. In recent works we showed that it is also a useful method to deal with 
quantum random spin systems with long range interactions, such as the Ising 
model with randon exchange interaction in the presence of a transverse and 
longitudinal magnetic fields [10] and the random Heisenberg model [11]. Most 
of the analytical and numerical methods to study these kind of systems rely on 
the so called replica-trick. Unfortunately, this clever technique becomes usually 
impractical within the glassy phase, where replica symmetry breaking occurs. 
As ED does not use replicas, this technical difficulty is not encountered, and 
the paramagnetic and glassy phases can be studied in the same way. Another 
appealing feature of ED is that it allows for the direct calculation of the 
dynamical response on the real frequency axis. In this way, it circumvents 
the uncertainties related to the analytical continuation procedures which are 
usually needed in quantum Monte Carlo simulations [14]. Finally, another 
important advantage of ED is the possibility of gaining insight on the nature 
of the low energy excitations. The price to pay is that only small systems 
are amenable to be treated within the available computer power. However, at 
least for the case of models with infinite range interactions, in most of the 
cases the relevant physical quantities were found to extrapolate smoothly to 
the thermodynamic limit and a consistent description of the different phases 
has been possible [10,11]. 

In this work we present results on the behavior of the anisotropic Heisen- 
berg model for a system of iV spins with S = 1/2 and random infinite-range 
exchange interactions. The Hamiltonian reads, 

1 N a 
H = -j= £ JiAStS! + -(S+S- + SrSf)], (1) 
viv i:j=1 * 

where i,j label the sites of the fully connected lattice and the interactions 
Jj j are normally distributed with variance J 2 that we set to unity. The pa- 
rameter a labels the degree of anisotropy interpolating between the classical 
Sherrington-Kirpatrick (SK) model for a = and the Heisenberg model for 
a — 1. This model has a rather close experimental realization in the above 
mentioned LiHo a .Y 1 _ :E F4 [4] compound in the limit of small magnetic ion con- 
centration x, which is dominated by random exchange interactions. The mag- 
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Fig. 1. The spectral function x'zi^) f° r the Heisenberg model (a = 1) in a system of 
N = 14 spins. The arrows indicate contributions from different kind of excitations 
(see text). 

netic interaction in this system is strongly anisotropic [8]. We calculate the 
different components of the local dynamical susceptibility, defined as 

1 m i N i 

xS» = ^E^ BClsf-^^sf iO, (2) 

1V1 M=l iV 1=1 LU U 

where m denotes the number of realizations of disorder and fi = x,y,z . 
The state |<&o"^) is the ground state of the Hamiltonian H^ m \ defined for the 
mth realization of Jij and it is calculated by recourse to Lanczos algorithm. 
The evaluation of the spectral functions x'Liv) = — 2Im[;\j t oc (y)] is achieved by 
using a continuous fraction representation of the dynamical response functions 
[5]. We perform averages over a number of realizations of disorder between 
M = 3000 for the largest systems (N = 15 spins) to M — 100000 for the 
smallest ones (N — 8 spins). 

Let us first summarize the main properties of the spectral function in the limit 
a — 1 (usual Heisenberg model) where the three components of the suscepti- 
bilities coincide [10,11]. This is shown in Fig.l. Four different contributions to 
the dynamical response can be distinguished: 

x» = ks{u) + x'LH + xL» + x'LM (3) 

The pure delta-function of the first term is a consequence of the SU(2) ro- 
tational invariance of the Hamiltonian in this limit. This response is due to 
a "soft mode" present in disorder realizations with total S ^ which have 
2S + 1 fold degenerate ground states. A close analysis of the finite size effects 
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on the average magnetization per site indicates that < S > /N — > as the size 
of the system evolves to the thermodynamic limit. Hence, this contribution to 
the spin response is expected to vanish as N increases. 

The analysis of the low frequency feature Xiow{ u ) as N grows reveals that this 
contribution of the spectral function gets sharper, evolving towards a func- 
tion ~ 8(u>). Its origin has been associated to slow coherent excitations of 
many-spin states which eventually become frozen. These states bear some re- 
semblance to the magnons of the model without disorder. In fact, by examina- 
tion of the structure of the ground-state in typical realizations of disorder one 
finds that these excitations are built-up on a ground state whose wave func- 
tion has large amplitudes on just a few configurations (out of the 2 N states) 
corresponding to unfrustrated sub-clusters, i.e. a sub-set of spins that are in a 
configuration that is compatible with the sign of the random bonds. Such 
configurations appear in pairs together with their time-reversal counterparts, 
having the same weight in the ground state wave function and some relative 
sign that depends on the total S of the ground state. The wave functions of 
excited states that contribute to Xiowi 00 ) & l so have a large weight on the un- 
frustrated clusters where the pairs of configurations appear with the opposite 
relative sign with respect to the ground state. The typical energy scale for 
these lowest energy excitations is 0(J/N). The high-energy part Xhighi 00 ) * s a 
small contribution with the form of a mild hump which is produced by exci- 
tations generated from the unbinding of single spins out of the unfrustrated 
clusters, the classical picture being a precession of individual spins around the 
effective quasi-static local field of the remaining frozen ones (of the unfrus- 
trated sub-cluster). It is remarkable that the three above mentioned features, 
which are the consequence of some kind of magnetic order are mounted on 
a broad and large background Xregi 00 ) that remains almost unaffected as the 
size of the system increases. This piece of the response has been identified as 
due to incoherent excitations. 

With this picture in mind, let us now carry out a similar analysis on the 
model with anisotropy. For a^O, SU(2) rotational invariance is broken and 
only the rotational invariance with respect to the z axis is preserved, so that 
Xx( u ) = Xy(u) 7^ Xzi 10 )- The behavior of these spectral functions for different 
values of the parameter a is shown in Fig. 2. Both spectral functions contain 
a piece of the form K8(u). Its origin is the degeneracy of ground states with 
large S z (or S x = S y components). In the case of x'x( u )i this piece carries 
a very low spectral weight and, in both cases, it is observed to decrease as 
the size of the system grows. Therefore, as in the pure Heisenberg case, this 
contribution is expected to vanish in the thermodynamic limit. Save from this 
feature, Xx(oo) does not show any clear indication of glassy order along this 
direction. The remaining part of the spectral function behaves as produced by 
an incoherent continuum, experimenting negligible changes when N is mod- 
ified, without developing low energy features. We identify it with the regular 
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Fig. 2. The spectral function x"( w ) (upper panel) and x"(u;) (lower panel) for 
a = 0.2, 0.4, 0.6 (bold line, dashed line, circles), in a system of N = 14 spins. 

part x'^regi 1 ^)- The spectral function corresponding to the response along the 
z direction also contains some portion of its weight distributed on a wide range 
of frequencies that is interpreted as caused by incoherent excitations. Inter- 
estingly, at least for strong enough anisotropy, Xx,reg( u ) an d x'^regi.^) seem to 
behave linearly in uo at low frequency, as predicted by a mean-field description 
which is exact for the SU(M) extension of the Heisenberg model in the limit 
of M— > oo [13]. A more detailed study of the origin of this behavior is left for 
future investigations. 

We now turn to analyze the contribution Xiow( u ) "which is clearly distinguished 
within the low frequency region of x'zi 00 )- m an analogous way to what hap- 
pens in the case of the isotropic Heisenberg model, the typical width of this 
feature is 0(aJ/N) and evolves towards a ~ 5(uj) as the system approaches 
the thermodynamic limit. Its origin is also due to excitations where the spin 
configurations conforming unfrustrated clusters, which in the ground state are 
frozen along the z direction, become canted with a finite component on the xy 
plane in a coherent fashion. Assuming that the functional form of the spec- 
tral function in the thermodynamic limit is x"(uj) = Q^i 00 ) + x'^regi 00 )^ the 
extrapolations for the weight accumulated in x'Lwi 00 ) as N — > oo provides an 
estimate of the Edwards Anderson parameter q. 

Recalling that the sum rule for the spectral functions is 
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Fig. 3. The integrated spectral weight corresponding to x"(u;) for a = 0.2,0.4,0.6 
(bold, dashed and dotted line) in a system of N = 14 spins. The arrows indicate 
the weight accumulated in the low frequency piece xiow{<^)- 



it is clear in Fig. 2 that an important transfer of spectral weight takes place 
from x'^regi 00 ) to <? as a approaches the classical limit a = 0. A more quanti- 
tative picture of this effect is shown in the logarithmic plot of the integrated 
weight I(uj) = ff^du'x'^uj') shown in Fig. 3. The weight accumulated in the 
low frequency peaks defines plateaus in I{oo) that are indicated with arrows 
in the figure. Extrapolations of the ordinates of these plateaus as the size 
of the system increases provides us a reliable estimate of the Edwards An- 
derson parameter q in the thermodynamic limit. Results are shown in the 
upper panel of Fig. 4. This plot can be reasonably fitted with the function 
q = 1/4^1 — (a/a c ) 2 , with a c ~ 1.02. Interestingly enough, this functional 
form suggests that, starting from the SK model, the fluctuations introduced 
via the interaction between the x, y components of the spins behave in a similar 
way to those driven by the temperature T in the pure classical model. In the 
latter case it has been found q — 1/4 oc T 2 [1], for T ~ 0. Instead, it seems that 
the effect of quantum fluctuations in the anisotropic Heisenberg model is differ- 
ent from those caused by adding a transverse field T to the classical SK model, 
where q = 1/4(1 — r/r c ) [15]. This kind of behavior should be relevant in the 
analysis of quantum vs classical annealing recently studied experimentaly in 
the LiHo x Y 1 _ a ,F 4 compound [7,8]. The energy per site as a function of the 
inverse of the system size is shown in the lower left panel of Fig. 4. The result 
of the extrapolations to the limit N — > oo are shown in the lower right panel. 
For comparison, we also show the estimates provided by the replica symmetric 
solution [16] and the one and two step replica symmetry breaking solutions 
[17] for the classical model, corresponding to e = —0.1995, —0.1913, —0.1909, 
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Fig. 4. Upper panel: The extrapolated values for the Edwards Anderson order pa- 
rameter as a function of the strength of anisotropy. Lower left panel: Energy per site 
as a function of the inverse of the lattice size. Lower right panel: Extrapolated en- 
ergy per site as a function of a. The square, diamond and triangle correspond to the 
replica symmetric [16], one and two step replica symmetry breaking [17] solutions. 

respectively [18]. It is worth to stress that, within that approximation, im- 
proving the level of symmetry breaking leads to higher values of the ground 
state energy per site. We have not explicitly investigated the limit a — 0. 
However, it is clear from Fig.4 that a simple monotonic continuation of the 
curve depicted by our numerical results leads at this limit to e ~ —0.183, i.e. 
to a value even larger than that predicted by the two-step replica symmetry 
breaking solution. As it is likely that higher levels of symmetry breaking give 
higher energies, we can say that our estimate provides an upper bound to the 
exact ground-state energy of the classical model. 

To conclude, we have investigated a model of a quantum spin-glass with 
long range anisotropic interactions which is a quantum generalization of the 
Sherrington-Kirkpatrick model. Our results show that both Xz( u ) an d x'xi 00 ) 
have a rich structure. We found that the z-component remains frozen for all 
values of the anisotropy and that the order parameter q decreases very fast 
for anisotropies a > 0.6, being remarkably weak in the pure Heisenberg limit. 

It is worth mentioning that for, large enough anisotropy, the low frequency 
edge of the regular part of Xxi 00 ) an d Xzi 00 ) seem to be consistent with linear 
behavior. This point is in contrast with our previous analysis for the SU(2) 
isotropic model [10,11], but it is in agreement with the results obtained in the 
large M limit [13]. 
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